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SUMMARY 


A wind tunnel test program was conducted on an eight foot 
diameter model rotor system to determine blade element air- 
loads characteristics in the unstalled and stalled flight 
regimes. The fully articulated model rotor system utilized 
three blades with a Vertol 23010-1.58 airfoil section, the 
blades being 1/7.5 scale models of the CH-47C rotor blades. 
Instrumentation was incorporated at the blade 75% radial sta- 
tion to measure pressure and skin friction distributions, 
surface streamline directions and local angle of attack. The 
test program was conducted in three phases; non-rotating, 
hover and forward flight at advance ratios of 0.15, 0.35 and 

0.60. 


Test data were analyzed with respect to providing insight to 
the mechanisms affecting blade stall, particularly retreating 
blade stall during forward flight conditions. From such data, 
an assessment was made as to the applicability of current 
theoretical analyses used for the prediction of blade element 
airloads in the stall regime. These analyses utilize empiri- 
cal corrections to account for portions of three dimensional, 
unsteady environment; namely, the oscillating airfoil motion 
and the swept wing effects. 

The test results have indicated that the blade stalling 
process is delayed in comparison with a static, two-dimensional 
airfoil. The stalling process appears to be a sequence of 
events that develop over a significant segment of the rotor 
azimuth, rather than a sharp and distinct breakdown of the 
flow field. These events can be identified as follows; 

1. Onset of separation-like disturbances within the 
boundary layer and a slight shift of the center of 
pressure . 

2. Beginning of a rapid rearward shift of the center 
of pressure and an increase in Cm* 

3. Collapse of the leading edge suction peak. 

4. Achievement of maximum Cn. 

5. Achievement of maximum Cm. 

6. Completion of the stalling process. 
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These observed stall events provide support to the hypothesis 
that retreating blade stall is a phenomenon of three major 
phases: (a) a dynamic overshoot of leading edge suction and 

Cr 4nax' a subsequent loss of leading edge suction and the 

shedding of circulation in some form, perhaps as a modified 
vortex, from the leading edge region, and (c) complete upper 
surface separation. 

Prior to the commencement of the stalling process, the test 
data indicate that the three-dimensional nature of the flow 
field does not significantly affect the blade element airloads. 
Rather, the forces and moments resemble those that have been 
observed on two-dimensional airfoils oscillating in pitch. 
Correlation of test data with the aforementioned theoretical 
analyses indicates that the current empirical theories that 
account for three-dimensional unsteady effects can adequately 
predict blade element airloado on the advancing blade and dur- 
ing the onset of stall, but not after complete blade stall has 
occurred. These shortcomings are attributed to uncertainties 
in the blade element angle of attack, the local dynamic 
pressure and in the basic behavior of unsteady three-dimen- 
sional separated flow fields. 

The analysis has drawn heavily on the skin friction gage data, 
none of which is presently in suitable form for publication. 
These data will be incorporated into a subsequent report 
available from Dr. W. J. McCroskey. 
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INTRODUCTION 


The application of quasistatic, two-dimensional airfoil data 
to describe the helicopter rotor blade environment has been 
shown by many researchers to be inadequate for prediction of 
rotor performance, rotor blade loads, and blade motion at the 
higher forward speeds and thrusts being experienced by the 
current generation of helicopters. With the advent of the 
more sophisticated digital and analog computers to aid in the 
solution of the rotor problem, effects characteristic of the 
three-dimensional dynamic environment (heretofore assumed 
small) can be included in rotor theoretical analyses. The 
nature of these effects, though, is not clearly understood. 

Tests previously conducted on an oscillating airfoil, (Refer- 
ence ( l) , and on an airfoil undergoing a rapid transient 
change of angle of attack, (Reference 2), have provided in- 
structive data on unsteady aerodynamic effects as well as con- 
siderable insight into the dynamic stall mechanism. Tests 
conducted on swept two-dimensional wings, (Reference 3), have 
shown that spanwise flow effects have a pronounced impact on 
the measured airloads. 

Attempts at empirically modeling these unsteady and three- 
dimensional effects have been made. Utilizing the principal 
of superposition, Harris, et al (Reference 4), took into 
account both the airload characteristics exhibited by a wing 
when swept and those characteristics associated with an air- 
foil undergoing pitching motion and generated an improved 
rotor theory (Figure 1) . These effects, when integrated over 
the rotor disc area, improved the correlation of rotor perfor- 
mance with test data. However, it was not established in that 
study that the predicted blade element airload characteristics 
were indeed representative of those in the rotor environment. 
Also, it was apparent that knowledge about the detailed be- 
havior of the flow in the boundary layer was needed. 

Therefore, the present investigation was undertaken to answer 
three questions: 

1. How do measured blade element airloads compare with cur- 
rent rotor blade load empirical methods which account for 
the three-dimensional unsteady environment? 

2. How realistically do the oscillating airfoil stall effects 
and swept -wing effects represent the actual mechanisms of 
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the rotor blade element stall? 

3. Have any basic factors been overlooked in the formulation 
of the rotor blade environment? 

To this end, an experimental program was designed to: 

1. Obtain measured blade element airload trends and represen- 
tative measured boundary layer flow characteristics in the 
rotor stall regime. 

2. Compare theory to test measurements of rotor blade air- 
loads and blade element angle of attack. 

3. Enhance the understanding of the mechanisms of blade stall. 

The experimental program consisted of testing an extensively 
instrumented three -bladed model rotor in both hover and for- 
ward flight regimes. Absolute pressure transducers, located 
at the 75-percent radial station, were used to measured blade 
element airloads. Skin friction gauges, located near the 
75-percent radial station on another blade, were used to 
measure boundary layer flow characteristics. Blade element 
angle of attack was also measured, using a differential pres- 
sure transducer. Blade root torsion and flap angle were also 
recorded. 
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TEST EQUIPMENT 


Wind Tunnel 

The test program was conducted in the Boeing V/STOL Wind 
Tunnel in Ridley Township, Pennsylvania. The tunnel, a 
closed-circuit, single-return type with a rectangular test 
section (20 x 20 feet), can be operated in three basic config- 
urations: solid, slotted, or open-^et. A maximum test sec- 

tion speed of approximately 260 knots is attainable in the 
solid-wall configuration. The circuit is cooled by an ex- 
change system with exhaust doors located upstream of the bell 
mouth and intake doors located downstream of the test section. 
Power to drive the tunnel fan is supplied by a 13,700-horse- 
power AC motor and a 1, 500 -horse power DC motor. 

Dynamic Rotor Test Stand 

The dynamic rotor test stand (DRTS) used during this model 
test was designed at the Vertol Division for use in testing 
single rotor configurations with blades up to 16 feet in dia- 
meter. The DRTS is sting-mounted for both hover and forward 
flight testing. The stand housing is pivoted such that the 
rotor shaft angle can be pitched forward 30 degrees and aft 
20 degrees by a remotely controlled actuator. Rotor system 
collective pitch can be remotely controlled in the range 0 to 
30 degrees, and the longitudinal and lateral cyclic pitch have 
maximum ranges of control of +10 to -6 degrees and +7 degrees 
respectively. A six-component strain gauge balance is incor- 
porated into the DRTS to measure steady and vibratory rotor 
loads. Power is supplied to the DRTS by either a 120-hcrse- 
power or a 4 00 -horse power motor with six gearbox reduction 
ratios available. A photograph of the model rotor, mounted on 
the DRTS in the Boeing V/STOL Wind Tunnel is shown in Figure 2. 

Rotor Blades 

The rotor blades used in the test program are geometrically 
scaled models of the production CH-47B blades. The blades 
were tested in a three-bladed rotor system configuration with 
constant chord and linear twist variation along the span. The 
airfoil section used was the Vertol 23010-1.58 shown in 
Figure 3. 
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AIRFOIL CONTOUR 




BLADE CONSTRUCTION 


FIGURE 3 VERTOL 23010-1.58 AIRFOIL 
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The blades are constructed in two basic subassemblies ; the 
spar and the aft section. Figure 3 provides a sketch of a 
typical blade cross-section. The spar is a "D" section of 
Scotch uniply fiberglass built up around a formed mandrel of 
foamed plastic which then becomes an integral nonstructural 
part of the blade. This subassembly is then cured under 
pressure at high temperature. The aft portion of the blade is 
comprised of balsa, which is joined to the aft section of the 
spar by a fiberglass wrapping around the entire assembly and 
bonded in a final assembly mold. 

Since chordwise pressure distributions were required in this 
program, certain structural modifications were incorporated 
to accommodate pressure sensing instrumentation. Depressions 
of .040 inch depth into the spar and .030 inch depth into the 
balsa aft section were formed at the 75-percent radial station 
to accommodate the extensive instrumentation to be located 
there. These modifications permitted the transducers to be 
mounted below the airfoil surface and hence, a smooth, proper- 
ly contoured airfoil surface could be maintained at this 
radial station. 


Basic data for the blades are as follows: 


Airfoil designation 

Chord 

Diameter 

Twist 

Root cutout 
Solidity 


- Vertol 23010-1.58 

- 3.336 inches 

- 8.0 feet 

- -10.8 degrees 

- . 192R 

- .067 


Pressure Transducers 

The transducers used for this test were Kulite LPL-125-5 semi- 
conductor diaphragm sensors, .125 inches in diameter with a 
pressure range of 0 to 25 psia, a frequency response in excess 
of 20,000 Hz, and a linear response to pressure throughout the 
range of interest. These transducers were extensively tested 
to determine their durability and acceleration and strain sen- 
sitivity. All transducers were calibrated at the same time in 
a pressure vessel. Figure 4 illustrates a calibration curve 
for a typical transducer showing the response of the trans- 
ducer through the test data reduction system. Even at very 
low pressure, the transducer accuracy is very good. 

Figure 5 provides a cutaway view of the transducer 
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FIGURE 5 


PRESSURE TRANSDUCER INSTALLATION 
IN THE ROTOR BLADE 
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installation, and a sketch of a pressure transducer. Table 
1 provides a listing of the transducer radial and chordwise 
positions on the blade and a listing of the data runs for 
which each transducer was operational. 

Initially, strain and acceleration effects were found to be 
present, so design changes to the transducer mounting and 
diaphragm were incorporated to minimize these effects. 

To alleviate the strain effects on the pressure sensing dia- 
phragm, the transducer, located in a depression on the blade 
at the 75 percent radial station, was mounted on an elasto- 
meric sandwich. To protect the diaphragm from foreign object 
damage, a screen was placed over the diaphragm. The area 
adjacent to the transducer was filled with a silastic rubber 
compound which covered the protective screen with the excep- 
tion of a .030 inch diameter hole through which air pressure 
could be sensed. This small part between the transducer and 
the air stream was desired so as to obtain a better approxima- 
tion to a point pressure measurement. To minimize the accel- 
eration effects on the transducer, a design change was made 
which moved the wire connections of the electronic load 
bridge on the transducer diaphragm from the center to the edge 
so their weight would not affect the diaphragm deflection. 
Wiring from the transducers to the slipring assembly on the 
rotor shaft was routed internally to maintain a smooth air- 
foil surface along the span and the depression around the 
transducer assembly was filled with an epoxy compound to 
maintain the airfoil contour. 

Skin Friction Gages 

Boundary layer measurements were taken on one of the blades 
using the heated-film skin-friction gage technique developed 
by Dr. W. J. McCroskey of the U. S. Army Air Mobility R&D 
Laboratory. This technique makes use of dual-f i lament gages 
that are photo-etched from thin nickel foil and bonded onto 
the surface of the blade in a manner similar to strain gages. 
In this test program, these gages were used to measure the 
skin friction coefficient and surface flow direction through- 
out the blade stall regime. This was to aid in understanding 
the nature of the boundary layer flow, to determine transition 
and separation, and to help identify the stall boundaries that 
were encountered. 

Heated-film skin-friction gages were mounted on both the 
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TABLE 1 

BLADE INSTRUMENTATION SUMMARY 


RT.AnP TNSTPnMFN'TA'T'TnM 


RUNS 

OPERATIONAL 

- 1 


PRESSURE RUNS 


r/R 

BLADE SERIAL 

NO. 


Absolute Pressure 


S/N 21 

S/N 22 

S/N 24 


Top A 

x/c = 

.041 

.754 



933-035 


Top B 

x/c = 

.058 

.750 



933-035 


Top C 

x/c = 

.102 

.750 



933-r J5 


Top D 

x/c = 

.158 

.750 



933-035 


Top E 



.750 



933-035 


Top F 


■H 

.750 



1 

cn 

Top G 



.750 



933-035 

w 

B 

Top H 

x/c = 

.617 

.750 



933-035 

o 

w 

Top I 

x/c = 

.702 

mm 



933-962 

*-4 

2 

Cj 

Top J 

x/c = 

.561 

.750 



933-962 

P 

w 

OS 

D 

to 

CO 

Tod K 

x/c = 

.970 

.750 



933-962 

Bottom A 

x/c = 

.033 

.754 



933-035 

s 

cu 

Bottom B 

x/c = 

.054 

.750 



933-035 


Bottom C 

x/c =_ 

.104 

.750 



933-035 


Bottom D 

x/c = 

.165 

.750 



933-035 


Bottom F 

x/c = 

.303 

.750 



933-035 


Bottom G 

x/c - 

.483 

.750 



933-962 

1 

Bottom I 

x/c 

.701 

.750 



933-962 

1 

Bottom K 

x/c * 

.955 

.750 



933-035 


Dil’rerentiai 

Pressure 

.78 



933-035 

3 5 

= =>2 

Skin Friction Mag 
x/c ■ 

.1 

.723 



Random 

2 S S 

H CO H 

Root Flap Angle 

.045 



000-t 035 

□ 2 Z 

S H i 

Root Torsion 

.2 



000-035 


11 






































upper and lower surfaces of one blade at approximately the 
75 percent radial station (Figure 6) . Refer to Table 2 for 
the measured radial and chordwise locations. 

In addition, a V-probe hot wire sensor was mounted approxi- 
mately 0.1 inches above the upper surface on another blade. 

A film gage was mounted directly on the blade surface beneath 
the hot wire sensor. This instrumentation was utilized to 
help distinguish between fluctuations and dynamic disturbances 
originating within the boundary layer and those in the in- 
viscid flow. These two sensors would, for example, help to 
distinguish between an abrupt flow change due to a blade- 
vortex interaction and the passage of a localized separation 
bubble lying completely within the boundary layer. 

The heated-film skin-friction gages operate on the basic 
principle that the local skin friction coefficient, Cf, and 
the local streamline direction at the surface of the blade 
are proportional to the rates of convective heat transfer 
from the heated films to the air stream. The technique is 
analogous to conventional hot wire anemometry, except that 
special circuitry was utilized to provide the accuracy re- 
quired. The complete system had sufficient dynamic response 
to measure and distinguish between periodic signals due to 
the various harmonics of blade motion and random fluctuations 
due to turbulence and separation in the boundary layer. 

Since the gages were permanently bonded to the surface of the 
rotor blade, it was not possible to perform individual cali- 
bration experiments on the sensors. Therefore, the calibra- 
tion constants for the magnitude of the skin friction were 
obtained with the aid of numerical laminar calculations of 
the boundary layer at low angles of attack in the nonrotating 
and hover portions of the test program. For these calcula- 
tions, the flow was assumed to behave in a conventional two- 
dimensional manner. The calibration constants for the flow 
angle determination were obtained by assuming that the mea- 
surements at low collective pitch in forward flight should 
follow approximately the theoretical inviscid flow angle in 
the absence of any local disturbances. A detailed discussion 
of the skin friction results obtained during this test will 
be published at a later date. 
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I ADDITIONAL 

I INSTRUMEN- j SKIN FRICTION TRANSDUCERS 


TABLE 2 

BLADE INSTRUMENTATION SUMMARY 


BLADE INSTRUMENTATION 
SKIN FRICTION RUNS 


r/R 


RUNS OPERATIONAL 


Top A Mag 


x/c = .1 ! .723 


S/N 21 S/N 22 S/N 23 


4 - 


036-055 


Top B Dir x/c = .2 j .729 J 


036-055 


Top B Mag x/c = .2 .729 


036-055 


Top C Mag 


x/c = .3 ; .736 


036-055 


Top D Dir 


x/c = .5 ; .750 i 


036-055 


Top D Mag 


x/c = . 5 ; 


Top E Mag 


x/c = .6 i 


750 
757 ! 


036-055 

036-055 


Top F Dir 


x/c = .85| 


774 




Top F Mag 


x/c = .85: .774 ! 


036-055 i 


036-055 


Bottom A Mag x/c = .10; .723 


036-055 i 


Bottom C Mag x/c « .3 ; .736 ; 


Bottom E 

Mag 

x/c = .6 j 

.757 

Bottom F 

Mag 

x/c = .85 

.774 

Bottom F 

Dir 

x/c = .85 

! .774 

Elevated 

Dir 

x/c = . 50 


Elevated 

Mag 

x/c = .50 

.75 


Surface Dir 


Surface Mag 


x/c 

x/c 


.50 

.50 


Differential Pressure 


Root Flap Angle 


s 


jrjRoot Torsion 

< 


.74 

.74 


78 


045 


Random 


036-055 


036-055 


036-055 


036-055 i 

1 

036-055 ' 


036-055 


036-055 


049-055 


049-055 
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Differential Pressure Transducer and 

Angle of Attack Measurement Technique 

Angle of attack measurements were made through a differential 
pressure transducer located at the leading edge of the airfoil 
at r/R = .78. The differential pressure was measured between 
an x/c = 0 and an x/c = .020 chordwise position on the lower 
surface of the airfoil. The method of determining the angle 
of attack from the differential pressure, proposed by the U.S. 
Army Air Mobility R&D Laboratory, was based on two observa- 
tions. 

A unique relationship exists between the pressure distribution, 
the aerodynamic angle of attack and the local Mach number at a 
point on an airfoil for two-dimensional, steady, inviscid flow. 
The flow in the plane of a rotor blade cross-section generally 
resembles the two-dimensional quasistatic flow past that air- 
foil section if the blade is unstalled and the spanwise gradi- 
ents are small. These idealized conditions are more nearly 
realized in practice near the leading edge of a blade than 
elsewhere, and at the leading edge, the pressure gradients are 
large. Therefore, the relationship between the differential 
static pressure at two neighboring points near the leading 
edge, the local dynamic pressure in the chordwise direction, 
and the local aerodynamic angle of attack in the plane of the 
cross-section of the blade can be expressed in some functional 
form 

p l- p 2 

ACp = IjpvScOS^A = T<«» <*» M > 

or 

a = ^ ( ACp , a , M) 

where pi and p2 are the static pressures at points 1 and 2, 
p is the density, V is the magnitude of the local total vel- 
ocity, A is the local sweep angle, a is the angle .between the 
chord line and the normal component of velocity, a is the time 
range of change of a , and M is the normal component of the 
local Mach number; that is, M based on V cos A. 

Differential pressure transducers for angle of attack measure- 
ment were used on both the pressure blade and the skin fric- 
tion blade. Calibration of the differential pressure trans- 
ducer was accomplished during the nonrotating portion of the 
test program. The airfoil was pitched through a series of 
angles of attack. During this pitch sweep, the differential 
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pressure and pressure distribution about the airfoil at .75R 
were recorded and the blade element angle at this station was 
measured using an inclinometer mounted on an airfoil template. 
Due to the inaccuracies involved in locating the leading edge 
point of angle of attack reference line with the template, a 
check of this calibration was made. 

Such a check could be made by comparing an^le of zero normal 
force points between present test results and standard V23010- 
1.58 airfoil characteristics. By integrating the measured 
pressures around the airfoil, a curve of normal force coeffi- 
cient versus measured blade element angle was made. Figure 10 
presents such a curve compared to test results of the same 
V2 3010-1. 58 airfoil at full scale Reynolds numbers. The 
difference in the angle of zero normal force, 1.7 degrees, is 
attributed to the difficulty in obtaining an accurate refer- 
ence line with the template. Since the full scale data is 
referenced to the NACA chord line, a correction of 1.7 degrees 
was made to the measured angle of attack data to adjust it to 
the standard NACA reference. All angle of attack measurements 
presented in this report have been corrected to the NACA 
reference line (Figure 3) . The correction to the measured 
angle of attack is included in the calibration of the differ- 
ential pressure transducers shown in Figure 11. The data 
obtained from the differential pressure transducer yielded 
angles of attack only up to where stall occurred on the air- 
foil ( a = 14 degrees). There, a rapid decrease in differen- 
tial pressure was noted (Figure 11) and the angle of attack 
obtained from the transducer was invalid until re-attachment 
of the flow occurred. 

Trends observed below stall were used to extrapolate the cali- 
bration to higher differential pressures; for when the trans- 
ducer was flown in forward flight, much higher differential 
pressures were observed than those obtained statically. 

These higher differential pressures were expected and can be 
associated with the dynamic extension of Cn observed on os- 
cillating airfoils. 

The effects of sweep, Mach number, and rate of change of angle 
of attack were neglected for the data analysis of this test. 
The sweep effect is illustrated in Figure 7 by results obtained 
on a NACA 0012 in the U.S. Army Air Mobility R&D Laboratory 
(7 x 10-foot) subsonic wind tunnel. The pressure pi was 
measured at the x/c = .022 location on the lower surface, and 
the pressure p2 was measured at the leading edge. These 
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results show that ACp is sensitive to angle of attack and is 
mildly dependent on sweep angle below the onset of stall. 

The Mach number effect is illustrated in Figure 8 by results 
obtained from measurements by Lizak (Reference 5) . The re- 
sults indicate a strong Mach number dependence, but because of 
the angle of attack and the Mach number environment of a rotor, 
the effect on a determined by ACp would be small (i.e., when 
Mach number is large at ip = 90, a is small; when Mach number 
is small at ^ = 270, a is large) . Mach number effects were, 
therefore, neglected. The rate of change of angle of attack 
( a) effect is illustrated by Figure 9 from results obtained 
by Liiva (Reference 1) during an oscillating airfoil test. 

The effect is similar to what would be expected based on poten- 
tial flow theory for oscillating airfoils (Reference 6). The 
effect of a on the angle of attack measurements is not well 
defined. However, the static and hover angles of attack de- 
termined by the differential pressure method do agree well 
with measured angles as shown in Figure 17. 

Blade Loads and Motion Gages 

Due to the limited number of available data recording channels, 
only root torsion and root flap angle data for the pressure 
blade were recorded during the test program. These data, 
along with the cyclic input and knowledge of the test condi- 
tion, were used to reconstruct the blade element environment. 

Data Acquisition System 

The Boeing-Vertol wind tunnel data system is capable of accep- 
ting 120 channels of model and tunnel data. It samples and 
conditions the desired data channels, converts the data to 
digital form, and then operates on the data with an IBM 1800 
data processor. The data are then stored on tape and then 
transferred to the IBM 360 for further processing. On-line 
data in final form are available through use of the IBM 1800 
data processor. A block diagram of the data flow is given in 
Figure 12. 

The signal conditioners accept data signals from the model and 
tunnel sensors and provide power adjustment, amplification, 
zero setting, and filtering for each channel. The conditioned 
data are then supplied to one of the two multiplexers, depend- 
ing upon the nature of the data. 
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The low-speed multiplexer is capable of accepting a maximum of 
48 data channels. It is used to sample data at rates from 
100 Hz to 10,000 Hz and is intended for use with steady-state 
data where harmonic content is not of primary concern. The 
remaining 64 channels are connected to a high-speed multi- 
plexer which is capable of sampling dynamic data on line with- 
out loss of its harmonic content. Both multiplexers contain 
analog-to-digital converts which transform the data to digital 
form. 

The IBM 1800 data processor is a combination of general pur- 
pose computer and special input/output features. The activity 
of the 1800 is under the control of a series of operating 
programs. These programs will, among other functions, cali- 
brate and check the calibration of the instrumentation, take 
wind off zeros, monitor and display selected data, and acquire 
and process the data for final output. 

For each data point, a record of approximately two seconds was 
made. From this record, the data were processed through the 
system described above and then input into the various manipu- 
lative and analytic programs, where the data were converted to 
engineering units, harmonically analyzed, corrected, and pre- 
sented in a convenient form (Figure 13) . 
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TEST PROGRAM 


Test Procedure 

Data during the test program were acquired during two separate 
test portions because of the extensive instrumentation used. 
During the first portion, the rotor was configured to include 
the pressure instrumented rotor blade along with two non- 
instrumented blades. Pressure data, as well as limited data 
on blade loads, blade motion, and skin friction, were obtained. 
During the second portion of the test program, the pressure 
blade was replaced by the skin friction blade and the skin 
friction data were acquired at the same test conditions as the 
pressure blade. Instrumentation for differential pressure and 
skin friction at x/c = .1 on the upper surface were included 
on both blades to evaluate the repeatability of the two sets 
of test runs. 

Each of the test portions was further subdivided into three 
phases: nonrotating testing, hover, and forward flight. The 

nonrotating phase of the test program, accomplished first, 
served to calibrate the differential pressure transducers as a 
function of angle of attack and to determine the static aero- 
dynamic characteristics of the airfoil section. These static 
characteristics were to be used in theory versus test correla- 
tion and for comparisons of pressure distributions at th same 
angle of attack in the nonrotating and rotating regime. j. 

During this nonrotating testing, the blade was positioned at 
right angles to the wind and supported at its tip and mid-span 
by a fixture. Runs in which the nonrotating pressure data 
were acquired are numbered from 930 to 960. 

The hover phase of the test program, conducted initially at a 
tip speed of 250 feet/second, acquired data through a collec- 
tive pitch sweep. Increments between each collective setting 
were approximately two degrees. Later in the test program, 
when skin friction data were being acquired, a hover run was 
also conducted at a tip speed of 500 feet/second. 

Testing during the forward flight phase consisted of setting a 
collective pitch and varying the rotor shaft angle. The shaft 
angle was increased until some blade or test stand limiting 
condition was reached. Testing was conducted at advance 
ratios of 0.15, 0.35, and 0.60 at tip speeds of 250, 400 and 
500 feet/second. Refer to Table 3 for a listing of test runs 
and the rotor operating conditions. 
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TABLE 3 

MATRIX OF TEST RUNS 






r 




V«p 

TEST 

* 

CUFF 

RUN NO 


TYPE OF RUN 



Pr 

Off 

2, 26 



o 

250 

Pr 

On 

3 , 27 


Hover Runs 

II 


SF 


50 


Collective 


500 

SF 


55 


Pitch Stfeep 

LO 

rH 

• 

250 

Pr 

Off 

16, 18 


Forward Fit 
Shaft Angle 
Sweep 

H 

[ 3 . 

Pr 

On 

17, 19 





Pr 

Off 

12,22,2^ 


Forward 


250 

Pr 

On 

13,23,25 


Flight 

li) 

cn 

• 


SF 


51 


Shaft 

II 

zi 


Pr 

Off 

28 , 29 , 3C 


Angle 

Sweep 


500 

Pr 

On 






SF 

i 


54 



■■ ■ i 


— 

Pr 

Off 

14, 20 


Forward 


250 

Pr 

On 

15, 21 


Flight 

o 

• 


SF 


52 


Shaft 

Angle 

SweeD 

ii 


Pr 

Off 

31, 33 


3. 

400 

Pr 

On 

32, 34 


ff W W 



SF 


53 



8 

II 


Pr 

Off 

933+960 


Nonrotating 

0 

Pr 

On 

933+960 


Angle of 
Attack 

ZL 


SF 


36+48 


Sweep 


* Pr = Pressure Measurement Portion of the Test 

SF = Skin Friction Measurement Portion of the Test 



It is well to note that prior to acquiring data at eacft test 
run, a tunnel "blowdown" was conducted. The blowdown allowed 
sufficient time for the temperatures throughout the model and 
tunnel to stabilize. This technique is standard tunnel opera- 
ting procedure and is done to prevent any instrumentation zero 
shift difficulties. 


Instrumentation problems 

Early in the testing of the pressure blade, it became apparent 
that the extensive efforts made prior to the test to eliminate 
strain effects on the pressure transducers had not been com- 
pletely successful. Strain effects, though not large, were 
significant during the nonrotating testing because of the low 
pressures being experienced by the transducers. A technique 
was then developed through which the strain effects could be 
eliminated. For each run that was made with the pressure 
blade, a corresponding run was made in which the dynamic pres- 
sures were kept from the transducers by a cuff which was 
wrapped around the blade. The cuff was hollow so that all the 
transducers could experience a common pressure. The cuff was 
vented at the trailing edge so that flexing of the cuff due to 
the change in airloads experienced around the azimuth would 
not affect the pressure inside the cuff. These "cuff-on" 
data, then, were data including the strain effects, the accel- 
eration effects, and any pressure existing in the cuff. The 
pressure in the cuff would not affect the integrated pressures 
though, for the pressure would be the same for all transducers. 
The cuff-on data, obtained at identical conditions to the 
cuff-off data, could then be subtracted from the cuff-off data 
to account for the effects of strain and acceleration. 

An additional problem arose because of the "cuff correction" 
technique. As mentioned previously, a blowdown of the tunnel 
was conducted to stabilize all temperatures and minimize any 
instrumentation problems associated with a temperature change. 
When the cuff was placed over the transducer, it kept the air 
from moving across the blade and a change in temperature 
occurred at the transducer due to the resistance heating in 
the transducer load bridge. This caused a steady shift to 
occur in the pressures measured and the steady shift was 
different for each transducer. Because this would affect the 
pressure distribution, it had to be removed. This shift in 
the steady pressure values was isolated by observing the trans- 
ducer steady readings prior to starting a given run, and then 
immediately following the startup. The steady readings 
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recorded by the transducers could then be attributed to blade 
strain effects, as the short time between startup and the 
first data point precluded much zero shift due to temperature. 
This was done for each "cuff on" run, and the results were 
averaged to produce a set of probable steady pressure readings 
associated with rotation. These steady pressures along with a 
correction for the static pressure in the tunnel then replaced 
the measured "cuff on" steady pressures which included the 
large temperature induced zero shift. 

A far more serious problem arose when the pressure blade was 
rotated for the first time. As the tip speed passed 300 ft./ 
sec. during an rpm sweep in hover, transducers on the aft por- 
tion of the airfoil began to fail. In all, six transducers 
were lost (x/c = .70, .86, .97 on top; .48. .70, .96 on bot- 
tom) . Detailed examination of the blade and the transducers 
indicated that the failure occurred in the wiring internal to 
the blade. No repairs were attempted due to the time required 
and the test was continued. No additional pressure instrumen- 
tation was lost during the remainder of the program. 

Data Reduction 

The digitized data that were acquired on line and stored on 
tape were input to the data reduction path shown in Figure 13 
after the test was complete. The data was converted to engi- 
neering units, scrutinized for bad data samples, harmonically 
analyzed (24 harmonics) , and averaged over a number of rotor 
cycles. The pressure data were then corrected and processed 
through the Y-54 data computation package (Appendix) . Modi- 
fied "cuff on" corrections, described in the previous section, 
were applied to all hover and forward flight pressure runs ex- 
cept those at high tip speed (500 ft. /sec., M = .35); at high 
tip speeds, the cuff could not be kept on the blade. Because 
the trade off of bending moment and centrifugal effects tend 
to keep the strain effects constant with increasing tip speed, 
and the measured pressures increase approximately as the square 
of tip speed, strain corrections at these higher tip speeds do 
not significantly impact final results. 

The Y-54 data reduction program (see the Appendix for a de- 
tailed discussion of its calculations) reconstituted the 
pressure harmonics for each transducer, integrated the results, 
calculated the characteristics of the blade element environ- 
ment, and plotted the results. 

The skin friction data were provided to the U. S. Army Air 
Mobility R&D Laboratory in harmonic form for further data re- 
duction. 
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Data Analysis 


A total of six absolute pressure transducers, three on the 
upper airfoil surface at x/c =0.70, 0.86, and 0.97, and three 
on the lower airfoil surface at x/c =0.48, 0.70, and 0.96, 
failed early in the test program. In order to obtain aero- 
dynamic coefficients based on the integration of the pressure 
data around the total blade element surface, an engineering 
approximation was applied to those portions of the airfoil 
where pressure data were lacking. After viewing many pressure 
distributions from previous test programs and assessing a 
theoretical prediction of the pressure distribution, it was 
decided to approximate the pressure distribution over the aft 
of the airfoil with a parabola. This parabolic distribution 
was applied to both the aft uppei and aft lower surfaces in 
the following manner: The parahia passed through the most 

aft measured pressure data poin .nd had an infinite slope 
(dCp/dx/c = oo ) at x/c = 1.0 whe: Cp v.-s defined as zero. The 

equations for these curves are as follows: 


For the upper surface, 

^(Vc>.6l7) ~ 

For the lower surface, 

CV(x/ c > .303) = 

The sketch below illustrates the shape of these curves. An 
evaluation of this approximation is presented during the dis- 
cussion of the nonrotating data. 


C ! ~ 

^(•617) « *3 83 


r f / - */c 

L f > 030$) \l .617 




Pressure data were reconstituted at 10-degree increments of 
azimuth from the harmonic data. Blade element aerodynamic co- 
efficients (Cn# Cc# Cm) were then computed. The computations 
were made two different ways. The first involved integra- 
tion of only actual measured test data on the forward portions 
of the blade element; i.e., zero pressure was assumed aft of 
x/c = .617 on the upper surface and aft of x/c = .303 on the 
lower surface. The second utilized the engineering ap- 
proximation described above to account for the absence of test 
data on the aft portions of the airfoil. The trapezoidal 
technique was employed for the integration of the measured 
test data. For all aerodynamic coefficients presented and 
discussed in the body of this report, the approximation has 
been used. A complete listing, in microfiche form, of the 
blade element aerodynamic coefficients per ten degrees of 
azimuth position, for each run and test point (with and with- 
out the approximation applied) , can be obtained by contacting 
br. W. J. McCroskey of the U. S. Army Air Mobility R&D Labora- 
tory at Moffett Field, California. 

Calculations based on classical rotor theory. Reference 7, 
using the measured blade flapping were made to determine the 
blade velocity environment and the local dynamic pressure 
(qLoc) • This quantity was then used to nondimens ionalize the 
measured pressures and the integrated blade forces and moments 
and produce the aerodynamic coefficients. 

Blade root torsional load data were utilized to construct the 
blade elastic twist at the 75-percent radial station. Because 
blade elastic twist occurs predominantly in the first elastic 
torsional mode, knowledge of the torsional load at the root is 
sufficient to approximate the blade elastic motion at the 75- 
percent radius. This motion, along with a knowledge of the 
rotor cyclic and collective pitch, was used to calculate the 
blade total feathering motion. 

A theoret ical angle of attack as a function of azimuth was 
computed from the expression 

<*« * 6 * 1W ' 


using the blade feathering motion and the blade element calcu- 
lated velocity environment. 
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A second measure of the angle of attack was made based on 
knowledge of the blade element static airfoil characteristics 
and the dynamic Cn. It was calculated as 


Cm„„ ~ C-K 


- 


’DVNAm\C 




(« =■ o) Static 


STATIC 


An angle of attack was also calculated based on the differen- 
tial pressure measured at the leading edge of the airfoil (see 
Instrumentation section for background discussion) . This 
angle of attack, designated °^ml» was then corrected for the 
velocity at the leading edge of the airfoil due to blade 
feathering motion. When this correction is applied, the angle 
of attack is designated as o^vnlc * 

Other blade element environment quantities were also calcu- 
lated. A complete listing of these quantities, as well as the 
details of their calculations, is enclosed in the Appendix. 
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TEST RESULTS AND DISCUSSION 


Nonrotating (Static) Testing 

Airfoil characteristics up to an angle of attack of 22 degrees 
were obtained during static testing of the blade instrumented 
with the absolute pressure transducers. The testing was per- 
formed at a Reynolds number based on chord of .36 x 10^ and a 
Mach number of .17. The angle of attack was referenced to the 
NASA reference line. 

The curves of Figure 14 provide a comparison of static airfoil 
characteristics for the present test results at a Reynolds 
number of .36 x 10 6 and two dimensional airfoil test data at a 
similar Reynolds number of .3 x 10^ and a full scale Reynolds 
number of 2.6 x 10^. The comparison of the two normal force 
coefficient curves at nearly equal Reynolds numbers indicate 
that a stall type different from that expected of the V23010- 
1.58 airfoil was evident in the pressure blade data. The in- 
stallation of pressure measuring instrumentation at the three- 
quarter radius could hiive changed the airfoil contour, result- 
ing in a change in measured stall characteristics. The grad- 
ual change in lift curve slope and the reduction in maximum 
Cn indicate that the pressure blade was stalling from the 
trailing edge. Leading edge stall is characteristic of static 
V23010-1.58 airfoil data. From the integrated pressure data, 
lift stall is indicated at a s 12 degrees and a maximum C n of 
1.0 is achieved. At an angle of attack of 14 degrees, charac- 
teristics of moment stall begin to appear. This difference in 
the angles for lift stall and moment stall again point out the 
possibility of a difference in airfoil contour, for the two- 
dimensional test data indicates that these stall angles nearly 
coincide for the V2 3010-1. 58 airfoil. 

Pressure distributions accompanying the integrated data look 
two-dimensional below stall (Figure 15) . However, once stall 
occurs, as evidenced by no further increase in Cn with a , the 
pressure distributions differ from those one might expect 
knowing the characteristics of the airfoil at full-scale 
Reynolds number (Figure 16) . There is a rise in suction on 
the aft portion of the airfoil, but the collapse of the 
suction peak near the leading edge of the upper surface of the 
airfoil never clearly occurs. Both the pressure distribution 
and the integrated results suggest that stall is occurring 
from the trailing edge at this low Reynolds number. This is 
contrary to previous experience with the V23010-1.58 airfoil 
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which indicates leading edge stall even at this low Reynolds 
number and again suggests the possibility of a change in air- 
foil contour on the blades. Tho significance of this change 
cannot be evaluated from available data. 

Because the nonrotating data were the only data obtained 
during this test where all transducers were operational, it 
was utilized to check the accuracy of the engineering approxi- 
mation for the missing data. Analysis of the nonrotating data 
was made with the approximation replacing pressures on the aft 
portion of the airfoil where the transducers had failed during 
the rotating tests. Figures 15 and 16 illustrate a comparison 
of pressure distributions with and without the engineering ap- 
proximation that was discussed in the Data Analysis section. 
Figure 17 illustrates the effect of the approximation on the 
blade element force and moment coefficients throughout the 
measured angle of attack range. The use of the approximation 
causes Cn to be underpredicted below angles of attack of 
5 degrees. Above 5 degrees, the Cn is overpredicted. At an 
angle of attack of 10 degrees, Cn is overpredicted by .04 and 
at higher angles, Cn is overpredicted by .06. The blade ele- 
ment pitching moment coefficient. Cm# is slightly overpredic- 
ted at angles of attack below 5 degrees and is underpredicted 
(more negative) by approximately .025 at higher angles of 
attack . 

While the engineering approximation causes a shifting of the 
values of these coefficients, the trends based on the approxi- 
mation essentially parallel those established by the measured 
data. The similarity of trends provides a great deal of con- 
fidence in the use of the approximation. All such coeffi- 
cients computed in the hover and forward flight regimes and 
used for discussion in this report make use of this approxi- 
mation. 


Hover Testing 

Hover testing was conducted during the test program to deter- 
mine whether the rotating environment affected the pressure 
distributions and boundary layer flow characteristics that 
were observed during the nonrGtating phase of the test program 
The testing was initially conducted at a tip speed of 250 ft/ 
sec where the 75-percent radius blade element velocity environ 
ment was the same as in thG nonrotating test. A collective 
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pitch sweep was conducted from 8.75 = 0 to e .75 = 26 degrees 
in increments of two degrees. Figure 18 illustrates some 
typical hover performance results. 

Two independent sets of hover runs were made at identical con- 
ditions on separate days to ascertain data repeatability. Not 
only did the data from the two hover runs agree (Figure 17) , 
but it also exhibited excellent correlation with the inte- 
grated nonrotating pressure data. This correlation also veri- 
fied that the angle of attack measurement technique was work- 
ing, for the correlation was obtained only through comparing 
the pressure distribution at the same indicated angle of 
attack. 

Figures 19 and 20 illustrate the correlation of the rotating 
and static pressure distributions. Only at the leading edge 
on the upper surface can any difference between the rotating 
and nonrotating data be seen. 

Figure 17 illustrates the integrated pres sure /aerodynamic 
force coefficient data correlation between the hover and non- 
rotating runs. The data correlation is good, and the large 
number of angle of attack points obtained during the hover 
run provide a very dense calibration of the airfoil charac- 
teristics . 


Forward Flight 

In the forward flight regime, data were obtained at advance 
ratios of .15, .35, and .60. This was done initially at a tip 
speed of 250 feet/second because of an early experience with 
pressure transducer failure at the higher tip speeds. As the 
test progressed, however, and no further failures were en- 
countered, the tip speeds were increased and data were ob- 
tained at a tip speed of 400 feet/second at the advance ratio 
of . 6 , and 500 feet/second at the advance ratio of .35. Data 
runs were made with shaft angle sweeps at a constant collec- 
tive pitch up into the rotor stall regime. 

Rotor Performance 

Rotor performance data, although presenting a doubly integra- 
ted picture of the rotor environment, 

2tt r 

Performance = (xYlJxJY 
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is symptomatic of what is happening at the local blade 
element. Rotor problems like the decrease in control sen- 
sitivity to collective pitch at high rotor thrusts, drop off 
of propulsive capability as advance ratio increases and rapid 
increases in power required as thrust increases are sympto- 
matic of local blade element stall. The rotor performance 
measurements from this test are discussed, therefore, as in- 
dicators of the type of blade element events that can be ex- 
pected on the rotor. 

From observation of the forward flight performance data ob- 
tained at advance ratios of y = .15, y = .35, and M = .6 at 
tip speeds of 250 feet/second, a decrease in the rotor lift 
curve slope (d(Cr/cr/da s ) signaling the onset of both lift and 
moment stall is apparent. This decrease in slope (Figure 21) 
affects the rotor sensitivity to control inputs at high gross 
weights and adversely affects the helicopter fJving qualities. 
The lower graph of Figure 21 illustrates that as the advance 
ratio increases, the onset of rotor lift curve slope deter- 
ioration occurs at progressively lower values of thrust. Once 
stall is present, a large rise in the rotor equivalent drag 
takes place (D e = P/V -X). Figures 22, 23, and 24 illustrate 
the lift-drag polar s at advance ratios of .15, .35, and .6 
obtained during this test. The rapid rise in equivalent drag 
can be clearly seen at a pp/o = .075, y = .35, and also at a 
Gp/a = .06, y = .6. The thrust coefficients where a large 
rise in equivalent drag is noted also correlate with the 
thrust coefficients where a decrease in the rotor lift curve 
slope occurs. 

The propulsive efficiency of the rotor dX/d(P/V) is also 
symptomatic of stall on the rotor disc (Figure 25). At an 
advance ratio of .15 and a Gp/o of .08, the propulsive effi- 
ciency of the rotor is nearly 100 percent, indicating little 
increase in drag at a constant thrust as propulsive force is 
increases At a p = .35, the propulsive efficiency has 
dropped to around 88 percent at a pp/o = .08, but the rotor 
can still propel effectively. However, at a p = .6, the rotor 
is unable to propel itself at all. There, the propulsive 
efficiency varies from nearly 100 percent to negative values 
but at no time is the rotor propulsive force sufficient to 
overcome the rotor blade's own drag. At the high advance 
ratios, the drag on the advancing and retreating blades, 
where little lift is produced, creates a large force that the 
rotor must overcome to propel itself. As the power and shaft 
angle are increased at constant thrust to try to obtain more 
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propulsion, stall occurs on the retreating blade portion of the 
rotor disc. The loss of lift there, due to the stall, occurs 
where the rotor would use it more effectively for propulsion. 
This stall is not only responsible for reduced propulsion, but 
it also contributes largely to the rotor power required. It is 
interesting to note that large down flapping, and in some in- 
stances droop stop pounding, occurred at ^ = 50 to 60 degrees 
for these high advance ratio high thrust test runs, a direct 
result of the stall on the rotor disc. 

Stall Events at an Advance Ratio of *35 . - A particularly in- 
formative set of data for discussion of the stall mechanisms 
observed during the forward flight segment of this test was ob- 
tained at an advance ratio of .35 and a tip speed of 500 ft/ 
sec. The run consisted of a series of test points obtained 
during a shaft angle sweep from a s = -27 degrees, Ct/o = 0 to 
a shaft angle of -1 degree, Ct / a = .15. Collective pitch was 
held constant at 0.75 = 13.5 degrees and the tunnel speed was 
105 knots. Figure 26 illustrates key performance trends ob- 
tained during the run. From these performance plots, the stall 
onset on the rotor disc occurs after a Ct /o = .08. 

Figure 27 illustrates the angle of attack obtained from the 
differential pressure transducers on both the pressure blade 
and the skin friction blade for three of the test points. This 
degree of correlation ensures that the pressure results and 
skin friction results can be discussed at the same blade ele- 
ment and rotor environment. These points consist of data ob- 
tained just prior to stall onset, in moderate stall, and in 
deep stall and occur at Ct/ct ' s of .073, .105, and .132, re- 
spectively. Figures 28, 29, and 30 show integrated pressure 
data for the three test points under discussion. All three 
points show trends characteristic of an articulated rotor with 
minimum Cn on the advancing side and maximum Cn on the retreat- 
ing side of the rotor disc. Note that Cnmax occurs around tp = 
240 degrees for all three points. 

The Cm versus plot associated with the test point at a Ct/o 
= .073 indicates Cm (on the retreating side of the rotor disc) 
values that, when compared with levels from static data (Figure 
14) , would be statically characteristic of moderate moment 
stall at an angle of attack of 17 degrees. However, the usual 
sharp drops in Cm associated with moment stall are not seen. 
Although the performance curves of Figure 23 show no evidence 
of stall, the maximum normal force on the retreating side of 
the blade exceeds the static C nmax at the model scale Reynolds 
number. This value of normal force, however, does not exceed 
the static Cn ma x value at the full scale Reynolds number. The 
state of the rotor, therefore, is not clear, but it is inter- 
preted that this point is unstalled in lift and showing signs 
of moment stall onset. 

The variation in static Cnmax with azimuth at model scale 
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Reynolds numbers (Figure 28 and throughout the remainder of the 
report) was obtained from Boeing Vertol Wind Tunnel Test 074- 
077 where the variation of Cnmax with Mach number was acquired. 

At a Ct/o of .105, Figure 29, the rotor is stalled on the re- 
treating side of the disc. This is verified by investigation 
of the rotor performance curves of Figure 23. Cm is achieving 
values beyond those characteristic of the static stall, and 
the static stall values of Cn at full-scale Reynolds numbers 
have been significantly exceeded (dynamic lift overshoot) . A 
rapid change in both pitching moment and lift is also apparent 
as the blade enters the stall region. 

At a Ct/o of .132, Figure 30, the stall is severe. Cm ap- 
proaches values of -.6 and Cn is experiencing its maximum at 
values of 2.8. The rotor performance shown in Figure 26 has 
deteriorated significantly from the prestall values. 

These three test points reconfirm that on a rotor in forward 
flight, static values of Cn m ax are greatly exceeded. They 
also illustrate that the static values of Cm are exceeded as 
well, and the higher the thrust, the more dynamic overshoot 
occurs. In stall, both Cnmax an< * ^mmax are bein? reached at 
nearly the same time. Thir disagrees with static results 
which indicate that the maximum Cn and Cm do not occur simul- 
taneously and could indicate a further change in stall type 
due to dynamic characteristics. 

Figure 31 shows the Cn m ax that was achieved for each test point 
of the subject run and shows the stagnation pressure coeffi- 
cients, Cp, that were measured at the same azimuth position. 

The measured values greater than 1.0 are not normally expected. 
They are, however, representative of the local q and indicate 
that it was somewhat higher than expected. The local veloci- 
ties in the azimuth region where the data were obtained and q 
theory was calculated are in the order of 250 feet/second. 

They consist of inplane velocities due to rotation of the 
blades as well as velocities due to downwash and blade flap- 
ping. No significant velocities due to strong tip vortices 
are near the blade in the azimuth region where the maximum 
values are obtained. 

Further review of the blade element airload characteristics can 
be achieved by looking at the integrated pressure results as 
shown in Figure 32. These plots of Cn and Cm versus azimuth 
indicate that the test point of Ct/o = .132 is more severely 
stalled than that of Ct/o = .105. The fact that the Cn and Cm 
values from the three test points collapse to form a single 
line in the range 10 degrees <_ <_ 210 degrees indicates that 

nearly identical pressure distributions exist. The dynamic 
extension of the lift and moment coefficients is again ob- 
vious . 
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FIGURE 32 
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wnen Cn in plotted versus the angle of attack obtained from 
the differential pressure (Figure 33) , hysteresis phenomena 
are observed. Similarly, the curve of Figure 34 represents 
the hysteresis effect in the stalled region. These curves 
exhibit the basic characteristics of the hysteresis loops seen 
during oscillating airfoil tests at higher Reynolds number 
(Reference 1) and are somewhat similar to the theoretically 
predicted unsteady 3-dimensional versus a loops illustrated 
in Figure 35 for the test point at a C^/a = .132. 

The C n versus a loops illustrate the dynamic extension of lift 
over the statically measured values and indicate similar air- 
load trends for the blade element at all test points of the 
run when the angle of attack is increasing. 

Stall Onset at an Advance Ratio of .35: A Series of Separate 

Events . - The integrated results discvsred so far present a 
general view of the trends observed in the measured data. They 
have shown the similarity of the airloads in the first and 
second quadrants of the rotor disc and have also indicated that 
in the third quadrant of the rotor disc, dynamic extensions of 
both lift and moment exceed their static 2-D values. A more 
detailed look at the data is now in order- especially in the 
region of stall onset from p = 200 degrees to 240 degrees. 

In the third quadrant of the rotor disc, as the blade sweep:: 
around the azimuth past ^ = 180 degrees and toward p = 270 
degrees, it experiences a sequence of several major identifi- 
able events that occur prior to complete blade stall. In the 
present experiment, these events occurred over an azimi th 
region of 40 degrees and climax with the achievement of maximum 
C n • Blade element pressure distributions, skin friction dis- 
tributions, and the C n versus C m plot for the test point at 
y 1 = .35, Vt = 500 feet/second, Cy/a = .132, a s ~ -3.0 degrees 
are used in the following discussion to illustrate the events. 

As the blade traverses the advancing blade azimuth region from 
^ = 20 degrees to f = 180 degrees, the measured surface stream- 
lines essentially follow the theoretical local sweep angles 
constructed from Ut and cos p . In the second quadrant, the 
blade element circulation per unit span inferred from the 
integrated measured pressures, r = 1/2 V^,oc C L increases 
rapidly. The measured pressure distributions remain similar to 
those measured statically as shown in Figure 36. As the blade 
passes p = 200 degrees, however, the similarity tc static data 
is gone and subtle changes begin to take place while the stall 
onset region is entered. The skin friction instrumentation 
indicates a separation-like phenomenon, characteristized by 
large-scale random flow oscillation in the boundary layer, and 
the surface streamlines on the upper surface suddenly start 
to turn radially outward with respect to the theoretical local 
sweep angle. The center of pressure, r.table near the quarter 
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FIGURE 35 THEORETICAL ROTOR SECTION AIRFOIL CHARACTERISTICS 
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chord since ip = 120 degrees, begins to shift aft. The pres- 
sure distribution, even though at a C n of 2.0 and past the 
sLatic stall value, is still exhibiting a very strong suction 
peak, and as the blade passes this azimuth, it experiences a 
slight increase in suction at mid-chord on the upper surface. 
This shows up on the Cn versus C m plot of Figure 38 as a slicht 
change in slope and the data depart from a trend they have h - 
continuing since passing ip = 120 degrees. 


At the azimuth of 210 degrees, the lift, lift coefficient, and 
circulation are still increasing. The aft drift of the center 
of pressure accelerates and the onset of a rapid rise in nose- 
down pitching motion occurs (Figure 38) . This event has been 
described as moment stall in Reference 8 


and is definitely a 


separate event occurring during the stalling process. A further 
increase in suction is occurring at the mid-chord region and 
the leading edge suction peak is still clearly present as shown 
in Figure 37. 


The suction peak at the leading edge starts to collapse at 
ip ~ 215 degrees, but the mid-chord suction continues to grow. 
This keeps both the pitching moment and the normal force co- 
efficients growing rapidly. The circulation continues to 
increase but at a slower rate than before, and r finally reaches 
a maximum value at |i = 225 degrees. 


Stall increases in severity as the blade progresses in azimuth. 
At a ^ of approximately 240 degrees, after a considerable over- 
shoot of the static lift, the maximum C n is achieved. The 
pressure coefficient at mid-chord on the upper surface has just 
passed its minimum value. The pitching moment coefficient is 
also near its maximum negative value, considerably higher than 
that achieved statically. A few more degrees further in azi- 
muth and the lift stall occurs because of the redaction in the 
suction at mid-chord. As these events are occurring, the local 
circulation and lift force are decreasing at large rates. 

The blade continues on its way, through a region of deep stall, 
past a ip of 270 degrees and on to a of 345 degrees. Varia- 
tions of the pressure distribution occur throughout this 
region, but suction at the mid- and fore-chord on the upper 
surface remains a dominant feature. At = 345 degrees, the 
blade passes in close proximity of a vortex shed by the pre- 
ceding blade and then passes through the wake of the hub. At 
this time, the pressure distribution quickly changes to one 
similar to two-dimensional static data with an established 
leading edge suction. Stabilization of the center of pressure 
and re-establishment of the boundary layer unstalled flow 
characteristics soon follow as the blade begins another revo- 
lution. 

When the pressure distributions are viewed as a continuous 
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occurrence, the events described during the stalling process 
can be associated with the shedding of circulation from the 
leading edge of the airfoil. The large adverse pressure gradi- 
ents experienced by the airfoil at its leading edge are far 
above those obtainable statically and probably will exist only 
as long as the angle of attack is increasing at a large rate. 
Once the rate of increase of a becomes small, the flowfield 
becomes unstable and breakdown of the suction peak occurs . 
Details as to the nature of the circulation shed from the air- 
foil leading edge are not known at present. However, flow 
visualization of the process indicates that the shed circula- 
tion does not exhibit the properties of a shed vortex. Until 
further definition for this event can be made, all reference 
to it will be made as shed circulation. The change in flow 
characteristics measured in the boundary layer at \p = 200 
degrees signals the beginning of breakdown and, as the suction 
peak collapses, circulation is shed from the leading edge of 
the airfoil. This shed circulation then appears to sweep back 
across the airfoil at a velocity significantly less than the 
local freestream velocity, causing the various stall events as 
described above to occur. This circulation shedding phenomenon 
has been observed for some time in 2 -dimensional airfoil tests 
of oscillating wings. It is described by Ham and Young in 
Reference 9 and further developed in Reference 10. In fact, 
data presented in Reference 2 arising from a rapid transient 
high angle of attack change from 0 to 30 degrees on a 2- 
dimensional wing are strikingly similar to the data being 
currently discussed and is illustrated in Figure 38 in a C n 
versus Cm format. 

The sequence of stall events made itself felt on t' e blade 
torsion loads and it is interesting to observe the blade tor- 
sion waveform shown in Figure 39 with the stall sequence in 
mind. In this figure, the blade elastic twist, A9 e , is line- 
arly proportional to the measured torsional moment at the blade 
root. At (p = 210 degrees, the torsion load begins a rapid 
nose-down elastic motion in response to the large Cm's being 
generated. The largest nose-down motion is reached at = 245 
degrees, coincident with the maximum C m . With the following 
reduction in C m , the blade begins a rapid pitch-up motion and 
stall flutter has begun. The fact that the blade torsion has 
followed the blade aerodynamic load at r/R = .75 indicates that 
a significant portion of the blade is experiencing this sequence 
at the same time as postulated by Tarzanin in Reference 11. 

A summary of the series of events leading to complete stall can 
be obtained through Figure 40. This map, constructed from both 
pressure and boundary layer measurements, consists of a cross- 
plot of events from all the test points within the test runs 
at = .35, Vt = 500 feet/second, and 875 = 13.5 degrees. The 
shaded portion of the figure indicates regions where the air- 
loads can be described as qualitatively two-dimensionally static 
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in nature and quantitatively as two-dimens ionally unsteady (like 
substall oscillating airfoils) . The unshaded region indicates 
regions of retreating blade stall where a classical descrip- 
tion of the flow does not suffice. Here the sequence of events 
that can be consistently identified as leading up to complete 
blade element stall is shown. These events are: 

1. Onset of separation-like disturbances within the boundary 
layer and a slight shift of the center of pressure. 

2. Beginning of a rapid rearward shift of the center of pres- 
sure and decrease in C n . 

3. Collapse of the leading edge suction peak. 

4. Achievement of maximum C n . 

5. Achievement of maximum C m . 

6. Completion of the stalling process. 

7. Re-establ.ishment of leading edge suction. 

8. Stabilization of the center of pressure. 

9. Re-establishment of normal boundary layer characteristics. 

The stall regions illustrated indicate that the rotor first 
experiences moment stall with the resultant onset of a rapid 
rise in blade torsional loads and then as the thrust is in- 
creased, the rotor experiences both moment and lift stall. The 
two separately identifiable stall events (Reference 8 ) exhibit 
themselves on both blade loads and rotor performance. 

The dat~. obtained during the p = .35, V<r = 250 feet/second run, 
exhibit trends similar to those at p = .35, V? = 500 feet/second. 

Because the measured pressures were a quarter as large, the 
lower tipspeed run required correction for strain effects and 
zero shift. Figure 41 illustrates the stall region for the low 
tipspeed run. Figure 42 illustrates the correlation of mea- 
sured C n and C m at a C<r/o = .097. The fact that the two run 
results correlate well indicates that the "cuff correction" 
technique is working well at this advance ratio. It also 
suggests that the stall events are independent of tipspeed 
within the range of tipspeeds tested. 

This series of stall events that was observed characterizes the 
stalling process at the advance ratio of .35. The events seem 
to support the basic contention of Ham (Reference 10) that 
retreating blade stall is a dynamic phenomenon of three major 
phases: 
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1. A dynamic overshoot of statically experienced leading edge 
suction. 

2. A subsequent loss of leading edge suction and the shedding 
of circulation from the leading edge region of the airfoil. 

3. Complete upper surface separation. 

This present experiment defines the boundaries of these phases 
more clearly and furnishes additional details of the flowfield 
development, especially during the onset of stall. 

Airloads at an Advance Ratio of .15 - A map of the stall char- 
acteristics of the rotor at an advance ratio of y = .15 and a 
tip speed of 250 feet/second is presented in Figure 43. Mo- 
ment stall characteristics indicated by the aft shift in center 
of pressure begin at rotor thrust coefficients of C*r/ a = .085. 
The blade loads obtained at thrust coefficients below this 
stall level are presented in Figure 44. The values of Cn at 
Qr/o = .08 are well below the static Cn max level. The Cm 
levels, although low, are indicating the onset of moment stall. 
Figure 45 shows the Cn and Cm trends obtained at the highest 
thrust tested (Ct/<? = .126), and Figure 46 shows the Cn versus 
Cm loop. The associated pressure distr ibtuions at Cr/a = .126 
are given in Figures 47 and 48. These distributions appear 
similar to the static airfoil pressure distributions obtained 
in the non-rotating portion of the test when compared at simi- 
lar angles of attack. A collapse of the leading edge suction 
and a sharp dropoff in lift with increasing angle of attack, 
usual indicators of lift stall, did not occur even at this 
high thrust level. An investigation of the normal force curve 
of Figure 45 indicates that the maximum Cn experienced in the 
stall region does not exceed the maximum static Cn value at 
full scale Reynolds numbers. The airload trends observed from 
these data indicate that the collapse of the leading edge suc- 
tion and the rapid changes in the pressure distributions asso- 
ciated with the shedding of circulation from the leading edge 
do not occur at this Reynolds number and advance ratio. 

Because of the low dynamic pressures present at this advance 
ratio, even the cuff correction technique was not sufficient 
for the removal of high frequency strain effects. For this 
reason, the harmonic content of the pressure data was trun- 
cated after uh Q 5th harmonic. 
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Airloads at an Advance Ratio of .6 - Illustration of the stall 
characteristics at this advance ratio was obtained at a tip 
speed of 400 feet/second and a Cm/ a = .09. Up to stall onset 
on the retreating side of the rotor disc, the pressure distri- 
butions of Figure 49 retain the appearance of unstalled condi- 
tions, with well established leading edge suction, even though 
the measured pressures as well as the integrated results 
(Figure 53.) are greater than those obtained statically. As 
the blade progresses further in azimuth. Figure 50, the col- 
lapse of the leading edge suction and the subsequent shedding 
of circulation from the leading edge takes place and the stall 
events similar to those experienced at an advance ratio of 
m = .35 occur. The rearv;ard drift of the shed circulation can 
be traced in the curves of Figure 50 by the movement of the 
higher pressure area on the upper surface aft with increasing 
azimuth angle. Figure 52 confirms through a plot of Cn versus 
C m that trends of the integrated airloads similar to those 
obtained at u = .35 exist at the advance ratio of .6. A simi- 
lar comparison of the stall events at these two advance ratios 
can also be seen in the stall map of Figures 53 and 40 at 
250 feet/second. At both these advance ratios, the airloads 
in stall are dominated by the shed circulation due to the 
collapse of the large leading edge suction obtained through 
high positive rates of change of angle of attack while the 
angle of attack itself is high. 

The magnitude of the local dynamic pressures on the retreating 
portion of the rotor disc may have affected the data at this 
high advance ratio. The calculated local dynamic pressures 
reach minimum values of .01 psi for the 250 feet/second tip 
speed run and .03 psi for the 400 feet/second tip speed run. 
These pressures are on the order of the accuracy of the trans- 
ducers. The variations in measured pressures and the high 
steady pressure coefficients of the low q region could, there- 
fore, be due to transducer error, strain effects that were not 
completely removed by the cuff correction technique, and/or 
oscillatory pressures inside the cuff. It is also possible 
that three dimensional effects are playing a larger role at 
this high advance ratio. Additional work on the pressures 
measured in the azimuth region between \ P = 240 and \ if = 310 
degrees is required to clarify the source of the pressure 
coefficients obtained. 
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Additional Comments on the Airloads in Stall 


A Theory-Test Comparison . - Whether or not airfoil data 
achieved from two-dimensional testing are applicable to the 
three-dimensional unsteady rotor environment when the rotor is 
stall d has been of serious concern to these predicting rotor 
performance and blade loads. To shed some light on this sub- 
ject, a theory-test correlation was undertaken with the data 
previously discussed at a y = .35, Cq}/o = .132, as = -3 
degrees, 6i s = -5 degrees, Vt = 500 feet/second, and V = 105 
knots . 

Two analyses were used? one analog and one digital. Both 
utilized the empirical techniques described by Tarzanin in 
Reference 1 1 for accounting for the three-dimensional unsteady 
aerodynamic environment of the rotor blade. The correlation 
was done in two ways. The first involved the correlation of 
the test data with the theory run at the rotor test condition. 
The second utilized a constant angle of attack schedule in the 
analysis and compared various representations of the aerody- 
namic theoretical model including quasistatic aerodynamic 
alone, quasistatic theory with sweep effects, quasistatic 
theory modified for unsteady aerodynamic effects, and the 
theory utilizing both the sweep and unsteady aerodynamic 
effects . 

The angle of attack comparison is shown in Fic vr i 54. This 
schedule of angle of attack was obtained from the angle 
measured through differential pressure transducer up to an 
angle of attack of 20 degrees, where the data from the AP 
transducer become unusable for a determination. Above this 
angle, some engineering judgment was used to construct a from 
the blade cyclic input, blade live twist, and theoretical 
blade flapping. The angle of attack was then broken down into 
e. Up, and Ut components. 

Figure 55A shows the results when the low Reynolds number air- 
foil data obtained for the pressure blade during the nonrota- 
ting section of the subject model test were used at the con- 
structed angle of attvick, along with the unsteady aero effect 
empiricism obtained from full-scale Reynolds number oscillating 
airfoil tests. Although the unsteady aerodynamic effects do 
account for a Cn increase over Cn max predicted when the 
static data are used exclusively, that max Cn is far short of 
the one measured at the test point under discussion. At the 
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Cn max azimuth region, the Reynolds number for the static and 
rotating data for this test point are the same. 

Figure 55B shows the results obtained when full-scale Reynolds 
number data are utilized. The quasistatic theory alone is 
still inadequate on the retreating blade and the addition of 
sweep effects makes little difference. However, the combined 
three-dimensional unsteady theory r esults are very encourag- 
ing, for prior to Cn max, the correlation is quite good, with 
theory predicting the C n versus \p up to C n max. This result 
indicates that the model rotor blade is experiencing airloads 
similar to those experienced on full-scale rotors. That is, 
even though the local Reynolds number on the model blades is 
an order-of-magnitude lower than the full-scale rotor blades, 
the unsteady aero environment of the blade in the region where 
Cn is above 1.0 is causing airloads similar to those charac- 
teristic of high Reynolds number data. 

This result could be explained by interactions between un- 
steady aerodynamic effects and Reynolds numbers effects. 
Erickson and Reding (Referenced) point out that static trail- 
ing edge stall can sometimes be converted to leading edge 
stall by dynamic effects. This could also explain why the 
data at u = .15 show stall characteristics similar to the 
present low Reynolds number static data, while the stall 
characteristics of v = .35 and .6, where d<*/dt gets increas- 
ingly larger, resemble dynamic stall with vortex shedding from 
the leading edge of the airfoil. 

After the occurrence of Cn max, the correlation deteriorates 
and the theory generally underpredicts the measured airloads. 
However, the brief rise in C n after stall, associated theo- 
retically (Reference 12) to re-establishment of a sufficiently 
high increase in a and da/dt because of blade stall flutter 
(Figure 30) is both predicted and measured. The fact that 
the general level of the measured Cn values obtained in the 
azimuth region from \ l> = 240 to f « 360 degrees are so large is 
contrary to most expectations, for those measured Cn's remain 
well above the static values of C n max, even for airfoils at 
full-scale Reynolds numbers. This lack of correlation can 
possibly be attributed to the lack of understanding of sweep 
and plunging effects on airfoils in stall, for it is in this 
region of the azimuth where plunging is large. Although much 
work has been done with airfoils oscillating in pitch, both 
in and out of stall, little work has been done with plunging 
airfoils install. What correlation with current theory has 
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been done indicates serious deficiencies when the airfoil is 
stalled. 

A comparison of measured normal force coefficients with theo- 
retically predicted normal force coefficients, run at the 
rotor test condition with no constraints on angle of attack, 
is given in Figure 56. In this figure, an improvement over 
the correlations shown in Fiaure 55 is observed in the seament 
of the rotor disc between the azimuths of 90 degrees and 210 
degrees. This improvement over the results obtained from the 
specified angle of attack method of correlation is directly 
attributed to the a effects on the differential pressure 
transducer incorporated for measuring the angle of attack. 

The measured a lags the trv'- a in th'j region of ip where a is 
increasing. The correlate -i obtained is shown in Figure 56. 

Aerodynamic Uncertainties - It should be finally mentioned 
that the airloads in the complete stall regime, where the 
theory-test correlation deteriorate <are sensitive to blade- 
element angles of attack and local velocities. The values of 
local velocity that were inferred from Cp max exceeded the 
theoretical values consistently in the stall regime and the 
values of angle of attack inferred from the differential 
pressure were unusable in the stall regime. Both plunging and 
pitching motion make up the classical angle of attack, and 
whether this is a unique parameter as far as the blade element 
airloads is concerned is questionable and has not been 
answered by this experiment. The results of this investiga- 
tion do indicate that local velocities were estimated satis- 
factorily prior to stall, and that the angle of attack is 
directly related to the blade element airloads. However, 
determination of Viocal# and the meaning of in stall remains 
a formidable challenge in experimental rotor aerodynamics re- 
search . 
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CONCLUSIONS 


The flowfield around a helicopter rotor blade in forward 
flight is generally acknowledged to be complex, unsteady, and 
three-dimensional. The results of this experimental investi- 
gation confirm some existing conclusions about the airloads of 
the flowfield, as well as permit certain new conclusions to be 
drawn about the nature of the flowfield and airloads experi- 
enced by the blades. 

1. On the advancing blade region of the rotor disc, as well 
as during the onset of retreating blade stall, the three- 
dimensional nature of the flowfield does not significantly 
affect the blade-element airloads. Rather, the normal 
forces and moments resemble those that have been observed 
on two-dimensional airfoils oscillating in pitch and 
undergoing transient changes in angle of attack. 

2. The rate of change of the local blade element angle of 
attack appears to be one of the primary factors that de- 
termines the departures from static theory on the ad- 
vancing blade and during the delay in the onset of stall 
on the retreating blade. It also appears that the pitch 
rates encountered in high-speed flight, where da/dt is 
high, tend to induce dynamic leading edge stall, regard- 
less of the static stall characteristics. 

3. Retreating blade stall is a series of separate events that 
build up over a significant segment of the rotor azimuth, 
rather than a single breakdown of the local flowfield. 
Previously observed distinctions between "lift stall" and 
"moment stall" occur because of the dynamic nature of this 
series of events as circulation is shed from the airfoil 
leading edge and passes over the airfoil during a finite 
time span. 

4. The dynamic nature of the stall events is probably the 
primary reason that the basic stall mechanisms appear to 
be increasingly dominated by circulation shed from 

the leading edge of the airfoil as advance ratio increases. 
However, this does not preclude strong unsteady effects at 
low values of y if rapid changes in angle of attack are 
encountered; i.e., maneuvering flight, blade-vortex inter- 
sections, or blade elastic excitations. 
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5 . 


Contemporary empirical theories that account for three- 
dimensional unsteady effects can adequately predict blade- 
element forces and moments in the advancing blade and 
during the onset of stall but not after the blade stalls 
completely. The shortcomings are related to uncertainties 
in the angle of attack, in local blade element dynamic 
pressures, as well as in the basic behavior of unsteady 
three-dimensional separated flows. 

The present experimental investigation has explored the mech- 
anisms of retreating blade stall at one spanwise location of 
the rotor blade a little more deeply but much remains to be 
done. A continued investigation of the data obtained during 
the subject test program, as well as a review of other exist- 
ing data, is proceeding. 
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DATA REDUCTION PROGRAM Y-54 
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DATA REDUCTION PROGRAM Y-54 

A computer program has been developed which reduces data ob- 
tained during the pressure blade portion of the subject test 
program to useful engineering parameters . Data from pressure 
transducers, blade motion sensors, and blade loads gages are 
input to the program in harmonic coefficient form, along with 
data characterizing the rotor flight condition. These data 
are then operated on to provide time histories of the rotor 
blade environment, including: calculation of the blade ele- 

ment velocities, pressure coefficients, angle of attack, blade 
position, and blade motion; integration of Lhe pressure data 
utilizing a trapezoidal technique, and plotting of the data. 


mtation 


The following calculations are made based on the harmonic in- 
put data and the test run condition data. To simplify presen- 
tation of these computations, the following ground rules shall 
be used: 


At any given test point, each piece of data can be represented 
by a series of harmonic coefficients. 

Let the steady term of the harmonic -eries be associated with 
a subscript (o) . 

Let the remaining (i) terms of the sin (s) and cos (c) coeffi- 
cients of the harmonic series have the (ic) and (is) sub- 
scripts . 

Therefore, the test data can be characterized by the data har- 
monic coefficients, 

J = <1© + ZI ( dls #,N cos 

L* I ' 7 


Let each data type be represented by an identification symbol. 

Root Flapping = £ 

lie ad Lag Angle - X 

Blade Root Torsional Moment = RT 

Pressure = fa 

k = subscript designating a particular x/c 
and y/c transducer location, k * 1 to 20 
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Additional data defining test conditions required in the cal- 
culations are input to the program along with the harmonic 
coefficients. 

I. Flapping Deflection and Motion Normal to Disc Plane 

^ s P© + ( &it SINty* <► COS ty) Degrees 

£ s - 7^5 Z ft i- (Pis COS iy - 0u SIN iy) -*■ Radians/Sec. 

it 

= Jgo .21 (ftO*(-Pu SIN iy cojsty)^ Radians/Sec . 2 

l-l 

II. Lead Lag Deflection and Motion in the Disc Plane 

F 3 To + .i Cfi* S.H cy + cos ly) ^ Degrees 

y - Tss- - 7S * £, (**■)( Tu ««* l f - Lc *ih ir) ~ Ft. /Sec. 

V = T5o 75,1 jt (««.)*(- fc* *iNiV- f ie cos l\|r) «. pt. /sec . 2 

III. Elade Element Feathering Motion-Elastic Twist 

- * 

s K [Rr 0 -t ,Z (*t 1s S| W *-y + cos iy)] Degrees 

. h 

©e s K [ Z (fti)(ftTi* COS y - Rt*< iy)j -w Degrees/Sec. 

* K [ Z^ft**) (*“ ^Ti.s*lNiy - Rric C0si y)J ^ Degrees/Sec. ^ 

K = Deg/in-lb derived from first 
modal properties 


Kt^ = root torsion moment harmonics 
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Total Feathering 

e = ©o + .7se r + e,c cos y + e is sin y + Degrees 

e = 7^ [ CQ-X e,c siNljr + ©,s cosy) + A? ] Radians/See. 
e = lib [(^') J (-©ic c °S1(c- e,s SINiy) + Se ] Radians/sec. 2 
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VI. Determination o* Pressures 


PT01 

PT02 




+ ( 7*0. is S'" ly + »* 0 , £e COS tv) -V PS I. 

n , 

+ £, \ SIN ty + H* ic cos ty) ~ p$,. 


PT 19 = 

P 'S« + ( r ‘i S'N ty 

t P '\c 

cos t y) ~ psi. 

DPOI = 

0P oi. + ( OP.,. s SIN ty 

+ OP„ 

le cos ty) ~ psi. 

PTEO 

^Trf.f ~ 



(dpoi) = 

n 

2., ( Jll )( 0P OH* cosi ^ “ 

°P., U 

siNty) - psi./sec. 

(dpoi) = 

n 2 

Z (#0 ( - DPo/ . S/N LTfr - 
IS ! <*S 

0P.,. t 

cos ty) - PS i /see? 


VII. Determination of Pressure Coefficients 


CP - 

K ” 

^ r ' K > ( l } a V * 

V288 ' > L0C 

K- 1,20 

DC P = 

(OPOI ) 

1 

{.iSs) f 3 Vloc 


(DCP) = 

(dpoi) 

1 

(sis) j° vj c 


(o'cp) = 

t • 

(dpoi ) 

i 




(m) f ^oc 
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VIII. Determination of Angle of Attack 


Measured: 


°C mj = A + b(dcf) + C (dcp) * + D (DCP) 3 
<Km, = (DC p) [ 8 + 2C(DCP) + 3D(DCP) Z ] 


= (OCP) [ 8 + 2C(DCP) + 30(DCP)* ] 
+ (DCp) 2 [ ZC + feD (DCP) ] 


A, 8 , 0,0 - Curve fit constants supplied 

IN PROGRAM INPUTS . 


m,c = ***, 


- 4«CSIN [IPMcXeos °<m,~)Ce) 1 

L Vloc. i 


Calculated: 


‘“‘calc = © + W ( Up /u T ) 

°*calc ~ © +- t ~ W p*^r) 

v LOC 

« m . = C N - Cmq 


^ n 0 - C ^ (® - O 

( dc *</a-). -- ^| kso . 
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Trapezoidal Integration 


APPENDIX 


For the trapezoidal integration technique, the airfoil is 
divided into various regions. 

a. for calculation of normal force and associated 
moments 


L.E. 


Upper Surface 



Lower Surface 


b. for calculation of chord force and associated 
moments 

C 



Front 1 Rear 


The integration takes place at each where data is reconsti- 
tuted and uses the pressure coefficients of each transducer. 

When the engineering approximation for missing data on the aft 
of the airfoil is utilized, the program makes a series of 
pressure coefficients as follows; 
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Program Instructions for Including Engineering Approximation 

Make following pressure coefficients at the following ( x/c , 
y/'c) coordinates: 


XZc xZc 

-.052 .303 

-.044 .5 

-.038 .6 

-.027 .7 

-.012 .85 

-.003 .95 

.00 1.0 

.00 1.0 

.003 .95 

.011 .85 

.026 .70 

.037 .617 

Include these pressure coefficients in the set for integration 

when the engineering approximation option is called. 

Program Instructions for Doing Trapezoidal Integratio n 

I. Upper Surface Program Instructions 

(1) Sort transducers in ascending x/c order (y/c>0): 
Include transducer at x/c = 0 made from the 
differential pressure. 

(2) Integrate (Cp) on upper surface from x/c = 0 to 
last available transducer. 

(3) Integrate [(Cp)(x/c)] on the upper surface to get 
moment due to upper surface. 

II. Lower Surface Program Instructions 

(1) Sort transducers in ascending x/c order (y/c < 0) : 
Include transducer at x/c = 0 made from differen- 
tial pressure. 

(2) Integrate (Cp) on lower surface from x/c * 0 to 
last available transducer. 


Lower Surface 
Cp V/c = 


Co / / - X /c 

P Vc *.303 )j-^rr 


LO U/£* 


Upper Surface 


Pv/c 


- Cp 


*/ c - .6/7 




- */c 


383 
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(3) Integrate [(Cp) (x/c)] on the lower surface to get 
moment due to lower surface. 

III. Front Surface Program Instructions 

(1) Sort transducer into ascending order of y/c from 
(y/c < 0) to (y/c > 0) for x/c £ Card input #0012. 

(2) Integrate (Cp) on front surface from -y/c to +y/c 
over all transducers. 

(3) Integrate [(Cp) (y/c)] on the front surface to get 
moment due to the front surface. 

IV. Rear Surface Program Instructions 

(1) Sort transducers into ascending order of y/c from 
(y/c < 0) to (y/c>0) for x/c > Card input #0012. 

(2) Integrate (Cp) on rear surface from -y/c tc +y/c 
over all transducers. 

(3) Integrate [(C p )(y/c)] on rear surface to get 
moment due to rear surface. 

Normal Force Coefficient Calculation 


Subtract upper surface integrated pressure coefficients from 
lower surface integrated pressure coefficients. 

C N = Cu LO wcr ““ ^ N vrpen 


Chord Force Coefficient Calculation 


Subtract right hand surface integrated pressure coefficients 
from left hand surface integrated pressure coefficients. 

” ^ c PftOMT ” sack 
Moments About Leading Edge (Nose up +) 

Subtract bottom surface integrated moment J Cp ( x /c) <K x /c) 
coefficients from upper surface integrated moment coefficients. 
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Add the above to the difference between the right hand inte- 
grated moment from the left hand surface integrated moment. 


m 


l.e. 


= C 




L.e. 


~ C 
upper 


L.e. 


+• C 


LOU)£H 


'L.E. 


- C 


FRONT 




L.E. 


OR CK 


Harmonic Integration Option 

A program option exists that provides an integration of the 
harmonic coefficients of pressure. 

I. Coordinates of the Pressure Transducers Used for the Har- 
monic Integration 


Upper Surface 

Lower 

Surface 

x/c 

y/c 

x/c 

y/c 

XX { 0 

-021 

XX { 0 

-.021 

.041 

.010 

.033 

-.039 

.058 

.022 

.054 

-.041 

.102 

.037 

.104 

-.045 

.158 

.045 

.165 

-.047 

.206 

.049 

.303 

-.052 

.303 

.050 

(.483 

-,044 

.513 

.043 

* 5.701 

-.027 

.617 

.037 

V .955 

-.003 

(.702 

.026 



* S .861 

.011 



{.970 

.003 




* These transducers will not be present for all cases. 

XX This pressure is constructed from a differential 
pressure . 


P(0,-.0Zl) r 

P (.033 , 0 3^> - 

DPO 1 

u 

9-° 

•PC .033, - o 3 <)\ ■ 

- dP 0 

It 

to 

P‘(.033,-.03<|) ls 

- dP is 

n c = 

PC.033, -.039). 

9 L C 

- dPlc 
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II. Determination of Pressures Normal to the Chord 

Integration is separated into top and bottom surfaces. 
Upper Surface ; 


U o = 1 -Pl d( x /c) 


o ^0 


“is = J J (* /c ) 


u u * f'r A J(Ve) 


Upper Suppace 

PNO X /c 
L = 


Lower Surf ace : 


L, 

^i.S 

Lie 


f % c J ( VC ) 


Lower Surface 
ano x /c 


Total Normal Pressure: 


z 

0 


L 0 ~ 

N tS 

= 


N i,c 

S 

^"tC “ U Lc 
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$ 


*> 


Til. Determination of Pressures Parallel to the Chord 

integration is separated into left and right surfaces. 


Left surface x/c £ Card input value 
Right surface x/c > Card input value 


Left Surface: 


Fo 

Pis 

Fie 


v / c 


v / ( 


MIN 

y /c *•** 

,/ n it <* ( Y /c) 

C MIM 


v/ c WAX 


y /c 


/ ’“A.c 


WIN 


Right Surface : 


Left* Surface 
7 \ and v/c 

L - l,Z 4 - 


Go 

&ls 

& lc 


V/c WAX 


% 


f 

Ml* 



d( Y /c) 


max 

y/ f *•*<. JC^e) 

V <= WIN 

r /c wax 


wf T *ic d( V /0 

"c WIN 


Total Chord P r essure : 


Right Surface 
and v /c 
i. * i ,Zh- 


* 



*V-. 



/ ‘ 





II 

c 

0 

0 

<3 

1 

u? 

c is ~ 

*is- <^s 

Cic * 

F cc * S'ic 
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IV. Determination of Moments Due to the Pressures: 

Integration is done in four parts: 

Upper surface 
Lower surface 
Left surface 
Right surface 

Upper Surface : 

= j'^.OVc) J(%) 

nu is . /l- Aif (*Ve)i(%) 
fiu lc = { -fkA*Vc)<X*/c) 

0 “V 

Lower Surface : 

- / *.(*Vc) <K%) 

"*•« = / )d( r /c) 

MLu = S n h Cv«)^( % ) 
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$ 


4 


i 





Left Surface : 

V/c '■wx 

MF. = f -Pa, ( y Vc) d (%) 

Vc miN 
V/c <"m»x 

"^s = J ^ a ( y Vc)d(V c ) 

Vc 

V/c 

mf 6c - J n tc ( V Vc)d(V/c) 

V/c W <N 

Right Surface : 

Yfc, 

- ( ^.( V */c)j(%) 

v/c 

'Vc *MX 

Mt * • ( r Ai 3 (“Vc) d(Vc) 

V/c miN 
V/ c 

M&i. c = J (y»/ c )j(y/ e ) 

V/- m« M 

Total Moraent Pressures 


M# = Mu 0 - ML 0 t MF 0 -Mfr 0 
Mis = Ma i4 - ML t , MF is - MG C9 
M ic * Wu 4C - MLu + - MG tc 
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V. Computation of Resultant Harmonic and Phase Angle 

«t = /("tO* + (Ricf 

<(>(.= Arc T„n 




/'* A 


' ' tv .’••• 



■ \ 
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